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Various  radio-burst  and  proton  spectral  parameters  are  considered  as 
correlation  variables  to  develop  an  improved  prediction  scheme  for  the  proton 
peak  flux  magnitudes.  Besides  the  time  integrated  ^radio  fluxes  at  Individual 
frequencies,  we^also  use,  a frequency  integrat ion f over  two  different  ranges. 

For  the  protons , .jiia  hii<fei  considered  besides^ the  peak  flux  of  protons  .^with  energ; 
greater  than  10  MeV)  their  average  energy  and  their  duration, ^atuL-alse  applied 
a solar  longitudinal  correction.  Our  results  lead  to  an  improved  prediction 
scheme  and  provide  a better  understanding  of  the  correlation  phenomena. 
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I.  INTRODUCTION 

The  concept  of  using  solar  flare  radio  spectra  as  advance 

predictors  of  significant  proton  events  was  proposed  several 
1-3 

years  ago  at  AFGL.  The  original  studies  defined  the  basic 

criterion  (U-shaped  radio  spectrum)  which  indicates  with  a high 

degree  of  probability  the  occurrence  of  a proton 

event.  Subsequent  studies  addressed  the  question  of  predicting 

the  magnitude  of  the  proton  spectrum,  as  measured  by  the  proton 

peak  flux  for  the  energy  channel:  E > 10  MeV,  using  time  integrated 

4 5 

radio  flux-density  ' , peak  flux  density,  mean  flux  density  and 

5 6 

duration  , and  flash-phase  time  integrated  flux-density  . We 

7-9 

have  shown  elsewhere  that  the  slope  of  the  proton  spectrum 
can  be  predicted  from  the  width  of  the  U of  the  radio  spectrum. 
During  that  study,  we  carried  out  an  intensive  analysis  of  the 
radio  and  proton  events  from  1966  to  1973,  and  it  became  clear 
that  there  were  broad  correlations  between  various  features  of 
the  radio  and  proton  spectra  which  might  perhaps,  lead  to  an 
improved  prediction  for  the  proton  magnitudes. 

In  the  present  study,  we  have  considered  a variety  of  radio 
and  proton  variables,  based  on  the  central  idea  that  the  energy 
(at  the  source)  emitted  in  the  radio  spectrum  can  be  expected  to 
be  well  correlated  with  the  energy  (at  the  source)  carried  away 
by  the  proton  spectrum.  The  variables  considered  are  described 
in  section  II,  and  the  preliminary  application  of  these  variables 
over  a limited  data  base  indicated  successive  improvement  in 
correlations.  Factors  which  led  to  considerable  improvement  were 
the  locational  correction  and  the  average  energy  parameter  for  the 
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proton  spectra,  and  the  restriction  of  time  integration  to  the 
duration  of  the  U-pattern  for  the  radio  spectra.  Various  combi- 
nations of  the  proton  variables  were  then  correlated  to  U-time 
integrated  radio  fluxes  at  five  frequencies,  and  also  to  the 
frequency  integrated  radio  fluxes  as  described  in  section  III, 
using  the  entire  data  base  1966-73.  A prediction  scheme  for 
the  proton  magnitudes,  based  on  these  results  is  given  in  section 
IV  and  comments  for  further  studies  and  implications  of  this  study 
are  given  in  section  V. 


II.  CORRELATION  VARIABLES 


In  previous  studies , the  commonly  used  proton  variable 
has  been  I1Q,  the  peak  proton  flux  in  the  channel  E > 10  MeV. 

A variety  of  radio  variables  has  been  tried  earlier:  the  peak 

5 4 5 

flux  at  a given  frequency,  the  time  integrated  flux  ' at  a 

given  frequency,  (including  the  post-burst  increase) , or  the 

flash  phase  time  integrated  flux0  at  a given  frequency,  typically 

ranging  from  606  to  8800  MHz. 

On  general  physical  grounds , one  can  expect  the  total  energy 
emitted  into  the  radio  spectrum,  and  the  total  energy  carried 
away  by  the  protons  at  the  source  to  be  well  correlated.  Thus 
instead  of  I^q,  the  peak-flux,  one  should  consider  additional 
factors  for  the  proton  correlation  variable,  that  would  effectively 
reflect  the  energy  content  of  the  proton  emission.  The  factors 
we  have  considered  are  discussed  below.  Similarly,  we  consider 
additional  properties  of  the  radio  emission,  to  better  reflect 
the  energy  content  of  the  radio  emission  at  the  source. 

Proton  Factors: 

1) . Average  Energy  Factor  p^.  If  the  integral  spectrum, 
observed  near  the  earth,  is 

I ( > E)  ~ E~B 

it  is  easily  shown10  that  the  total  energy  carried  by  protons 

with  energy  > Eq  is  Eq  p^I(  > Eq) , 

where 


is  the  average  energy  factor.  Two  proton  events,  with  the  same 
peak  flux  I.Q  = I ( > 10  MeV),  would  have  different  energy  contents 
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depending  on  the  spectral  slope  B.  We  introduce  the  energy  factor 
by  considering  I^q  rather  than  I^q  as  our  correlate  variable. 

2)  Duration  t . The  total  energy  in  the  proton  emission 
clearly  depends  on  the  duration  of  the  event.  Events  with  the 
same  I^q  but  different  durations  will  have  differing  energy 
contents,  roughly  proportional  to  the  time  of  duration  t . We 
introduce  the  duration  factor  by  considering  I^qt  rather  than  1^0' 

3)  Solar  Location  A.  Two  events,  which  give  the  same  proton 
emission  at  the  solar  source  will  nevertheless  be  perceived  quite 
differently  at  the  earth  due  to  the  propagational  differences 
according  to  their  solar  longitudes.  Such  considerations  have 
been  used  in  other  related  studies and  we  apply  a solar  longi- 
tudxnal  correction  e , where  A in  radians  is  the  magnitude  of 

the  angular  distance  from  the  standard  reference  longitude:  57°W. 

. . 3A 

We  introduce  the  longitudinal  factor  by  considering  I^e  rather 

than 

In  the  pilot  phase,  we  examined  the  effect  of  introducing 
one  factor  at  a time,  for  a limited  data  base.  Typical  results 
for  correlation  with  the  radio  time  integrated  flux  at  8800  MHz 


were : 


variable : 

correlation: 

coefficient 


■■■io 

0.56 


I10pl 

0.61 


I10T 

0.62 


Z10P1T 

0.64 


The  data  sample  was  restricted  to  a set  of  neighboring  West-longi- 
tudes, so  as  to  reduce  the  variations  due  to  the  longitudinal 
factor.  It  was  clear  that  introducing  or  r produces  some 
improvement,  and  using  both  factors  provides  a greater  improve- 
ment. A similar  pilot  study  for  the  longitudinal  factor  showed 


that  the  improvement  in  correlations  was  quite  dramatic 
Radio  Factors: 
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1)  Range  of  Time  Integration.  In  earlier  studies,  the 
full  range  of  time  for  a given  event,  including  the  post  burst 
increase,  if  any,  was  used  to  obtain  the  time  integrated 

1(a))  = J*dt  I(w,t).  Some  events,  however,  may  be  compound 

events,  where  one  set  of  peaks  (at  different  frequencies)  is  at 
the  time  of  the  U-shaped  peak  flux  profile,  while  another  set  of 
peaks  occurs  much  later  (or  earlier)  and  is  not  part  of  the  U- 
phenomenon.  In  such  cases,  one  should  not  include  the  contribu- 
tions due  to  the  non-U  portions.  We  have  carefully  analyzed  our 
data  base  and  calculated  the  I(w)  by  restricting  the  range  of 
time  integrations  to  the  bursts  in  the  vicinity  of  the  time  of 
the  U-event.  We  found  this  to  be  especially  important  at  the 
two  lower  frequencies:  606  and  1415  MHz,  as  can  be  seen  from 
the  improvement  in  correlations  with  I^g  indicated  below. 

Range  of  Time  Integration 
all-time  U-only 

606  MHz  0.323  0.536 

1415  MHz  0.418  0.601 

2)  Integration  over  Frequencies.  In  addition  to  the  time 

integrated  I(w)  at  various  frequencies,  we  also  considered  the 
frequency  integrated  fluxes  Jl(w)  dw  over  the  following  ranges: 
(i)  606-8800  MHz,  since  these  are  easily  available  and  give  a 
standard  representation,  (ii)  to  , where  is  the  fre- 

quency of  the  observed  minimum  and  is  the  frequency  of  the 
peak  of  the  high  frequency  branch  of  the  U.  These  are  measures 
of  the  energy  of  radio  emission  over  the  corresponding  frequency 


ranges . 
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III.  CORRELATIONS 

The  correlations  between  the  various  variables  proposed 

in  the  preceding  section  are  described  here . The  data  base 

used  is  the  same  set  of  24  radio  U-Proton  events  we  studied 
. . 9 

earlier  in  another  context  , ranging  over  the  period  1966-73. 

In  view  of  the  radio  absorption  problems  at  the  limb,  3 events 

which  were  > 85°E  or  > 85°W  were  omitted.  Also,  one  event  was 

too  close  to  sunset  and  the  time  information  necessary  for  radio 

data  was  not  fully  available.  One  additional  event  did  not  have 

observations  at  some  of  the  intermediate  frequencies.  This  left 

19  events,  as  described  in  Table  I,  which  provides  the  date,  the 

location  and  type;  I1Q,  A,  p-^,  t for  the  protons  and  I (w)  for 

the  radio  data  at  five  frequencies,  from  606  to  8800  MHz. 

The  actual  variables  to  be  used  for  correlations  are  presented 
• 3A  3A 

in  Table  II:  I^Q,  I10  p1'  I10e  ' "^lO6  pl'  for  the  Proton  data 
and  I ( u) ) for  5 frequencies  and  fl  ( oj ) du>  over  two  frequency 
ranges  for  the  radio  data.  We  have  omitted  the  proton  duration 
study,  since  it  is  an  after  the  fact  observation  and  cannot  be 
prejudged  in  terms  of  making  a real  time  prediction  based  on  radio 

data.  This  objection  does  not  apply  to  the  determination  of 
. . 7-9 

since  it  can  be  predicted  from  the  width  (oj^/u^")  of  the  radio 
U-spectrum. 

The  correlation  coefficients  for  the  7 radio  vs.  the  4 proton 
variables  are  described  in  Table  III.  From  these  results,  we  can 
draw  the  following  conclusions. 

(1)  The  location  factor  e3A  provides  an  improvement  in 
almost  every  case.  I1()e3A  is  better  correlated  than  I1Q,  and 
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3A 

I10e  p1  is  better  correlated  than  I1Q  , for  every  radio 
variable  with  the  sole  exception  of  the  individual  frequency 
4995  MHz. 

(2)  The  energy  factor  p^,  as  in  I^g  P^^  , also  provides  an 

improvement  over  the  direct  I1Q  in  most  cases.  However,  the 

”*A  3A 

choice  between  I^e'  and  I^ge  P1  depends  on  which  radio  variable 
is  being  used  as  the  predictor. 

(3)  Amongst  the  various  radio  variables,  the  individual 

frequencies  show  a moderate  variation;  with  2695  MHz  the  best  for 
3A  3A 

I1Qe  and  I1Qe  p^  f and  4995  MHz  the  best  for  I1Q.  The  lower 
frequencies  606  and  1415  MHz  are  less  reliable  than  the  higher 
ones,  in  every  case. 

(4)  It  is  better  however,  to  use  a variable  integrated  over 

the  frequencies,  such  as  dw  over  606  to  8800,  since  it  will 

be  less  susceptible  to  the  variations  occuring  at  individual 
frequencies.  Also,  it  will  be  possible  to  obtain  such  an  inte- 
grated variable,  even  when  the  individual  frequencies  available 

are  not  the  ones  used  in  this  analysis.  In  fact,  the  correlations 
3A  3A 

f°r  I^ge  (r=0.806)  and  I^e  p^(r=0.798)  vs.  the  integrated 
(606-880C  MHz)  radio  flux  are  better  than  those  for  the  individual 
frequencies . 

The  best  fit  straight  line,  using  y = log  (I1Qe  A)  and 
x = log  di*>  ) from  606  to  8800  MHz  is  given  by 


a-  0.7110 


(1) 
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Taking  the  exponential  of  eq.  (1)  we  find  the  predicted 

Il0e3A  = (0.1146).  (/ x U)du)  )1  • 770  .(5.14)±1  (2) 

The  best  fit  straight  line,  using  y = log  (IlQe3A  p1)  and 
x = log  ( J I ( ) da)  ) from  606  to  8800  MHz  is  given  by 

Y=ax  + b±o  (3) 

a = 1.6548  a = 0.6831 

b = -0.1193 

Taking  the  exponential  of  eq.  (3) , we  find  the  predicted 

I10e3A  P1  = (0,7598).(/i  (w)  dw  )X  ' 6 54  8 . (4 . 82r 1 ( 4 ) 

The  data  points  and  the  best  fit  straight  lines  for  I1Qe3A 
3A  f— 

and  l^ge  vs  J I (ui)  do)  from  606  to  8800  MHz  are  represented 

graphically  in  Figures  1 and  2 respectively. 


VARIABLES  FOR  CORRELATIONS 
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TABLE  III 

CORRELATION  COEFFICIENTS 


Proton  Variable 


Radio 

Variable 

I10 

x 3A 

I10e 

3A 

'lO61  C1 

I10  1 

A.  I (u)  at: 

8800  MHz 

0.718 

0.758 

0.773 

0.727 

4995  MHz 

0.781 

0.776 

0.780 

0.780 

2695  MHz 

0.759 

0.777 

0.794 

0.771 

1415  MHz 

0.601 

0.749 

0.711 

0.546 

606  MHz 

0.536 

0.707 

0.662 

0.474 

B.  fi  (u)  du  over: 

606-8800 

0.753 

0.806 

0.798 

0.736 

u2  - “3 

0.502 

0.631 

0.712 

0.578 

— ■ 


TABLE  IV 


Px  vs  6 


6 

* 

P1 

6 

P1 

5.0 

1.25 

1.13 

7.5 

4.0 

1.33 

1.0 

8.0 

3.0 

1.50 

0.9 

8.5 

2.5 

1.67 

0.75 

9.0 

2.0 

2.00 

0.6 

10 

1.6 

2.25 

0.5 

11 

1.6 

2.67 

1.5 

3.00 

1.4 

3.50 

1.3 

4.33 

1.2 

6.00 

1.166 

7.00 

P1  6-1 

in  this  column. 

p^+oo  for  6 + 1. 

We  do 

formula  for  6 <1.166. 

**  p^,  as  a physical  energy  factor  must  remain  finite.  The  slope 
6 < 1 is  generally  valid  only  over  a limited  energy  range,  and 
then  the  true  P^  will  not  become  infinite.  We  have  assigned 
these  values  of  p^  for  6< 1.166. 
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IV.  PREDICTION  SCHEME 

From  the  correlations  described  in  the  preceding  section, 

we  can  develop  a real  time  prediction  scheme  for  obtaining  the 

proton  flux  magnitude  I1Q  using  radio  variables  and  the  location 

of  the  flare  as  input  variables.  From  Table  III,  we  observe 

that  either  the  individual  I(w)  at  2695  or  4995  or  8800  MHz,  or 

the  integrated  results  J*I(u))du)  from  606  to  8800  MHz  can  be  used 

as  the  input  radio  variables.  Out  of  the  four  predicted  quantities 

3A  3A 

I10'  J10  Pl'  I10e  ' Iioe  p l the  latter  two  generally  have  the 
better  correlations.  We  will  use  the  integrated  radio  variable 

606-8800  MHz,  since  it  provides  the  highest  correlations.  The 

. . 3A  3A 

’predicted’  quantities  I^q  e or  I^q  e are  given  respectively 

by  eqs . (2)  or  (4).  To  obtain  l^g,  one  must  multiply  by  e-3  for 
-3A  -1 

eq.  (2),  or  by  e p.^  , for  eq . (4). 

The  location  factor  A can  be  obtained  directly  from  the 

flare  position  data,  in  real  time,  and  that  correlation  is  easily 

-3A 

applied.  Thus  eq.  (2) , followed  by  the  e correction,  provides 
a reasonable  scheme  for  the  prediction  of  I^g. 

The  energy  factor 


P1  = 


(B>1) 


cannot  be  obtained  a priori.  However,  our  earlier  studies3  have 
shown  that  the  relation 

B . («) 

(logjjto.A;))1,01 

provides  a good  estimate  for  the  proton  spectral  slope  and, 
consequently  for  P.^  through  eq.  (5).  That  equation  is  valid 
only  for8>  1.  Should  one  obtain  a wide  U,  eq.  (6)  will  predict 
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3<1  and  then  one  must  supplement  eq  (5)  with  the  following 
empirical  relationship.  In  fact,  such  a relationship  was  used 
in  studying  the  correlations  described  in  the  previous  sections, 
where  5 of  the  events  in  Table  I had  directly  measured  (3<1.  The 
values  for  these  events  were  approximately  set  according  to 
Fig.  3 or  Table  I***-  which  departs  from  eq.  (5)  for  0< 7/6  or  p1  > 7. 
Thus  eq  (4) , followed  by  the  e correction  determined  from 

the  locational  and  radio  data,  using  the  procedure  described  above, 
provides  another  way  to  predict  I1Q. 

We  illustrate  both  of  these  procedures  by  a concrete  example: 
On  30  April  1976  , a 2B  solar  flare  occurred,  reaching  maximum 
phase  at  2110  UT.  From  the  real  time  radio  and  optical  data,  we 
find: 

Position:  S09W47  , A * 0.175  Rad 

x =/l  (w)  du>  from  606  to  8800  MHz:  64.0 
<*>2  = 2695  MHz,  u>3=  11641  MHz , u)3/u>2"  = 4.32. 

(i)  Applying  eq.  (2) , we  predict 

I10e3A  = (0.1146)  (64. 0)1*770  = 180.4 
. - 3 a 

and  applying  the  e correction  with  A = 0.175  provides 

-2  -1  -1 
I1Q  = 106.9  cm  sec  ster 

(ii)  Applying  eq.  (4) , we  predict 

I10e3A  P1  = (°-7598)  (64. 0)1-6548  «=  740.6 

To  determine  the  p^  correction,  we  use  eq.  (6)  which  gives  the 
prediction^  * 1.87,  and  from  Table  IV  or  Figure  3,  P^55 


2.15 
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-3A  -1 

Applying  the  e correction  to  the  above  equation,  we 

obtain  the  predicted 

—2  —1  —l 

I10  = 204.5  cm  sec  ster 

The  actual,  observed  I1Q  for  this  event  was  370.0  cm”2  sec"1  ster"1. 
Thus  both  methods  provide  a reasonable  prediction  in  this  case, 
well  within  the  standard  deviation  factors  of  (5.14),  and  (4.82), 
respectively. 
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V.  CONCLUDING  REMARKS 

In  this  study  correlating  radio  and  proton  variables,  with 
the  aim  of  predicting  the  proton  flux  magnitude  I^Q,  we  have 
shown  that: 

(1)  the  time  integration  over  the  radio  data  should  be 
restricted  to  the  time  of  the  U-only; 

(2)  it  is  better  to  use  an  integration  over  a frequency 
range  such  as  606-8800  MHz  rather  than  use  I (co)  at  individual 
frequencies ; 

3A 

(3)  the  locational  factor  e improves  the  correlations; 

(4)  the  proton  average  energy  factor  p sometimes  improves 
the  correlations. 

Other  factors  that  were  considered,  but  not  fully  explored, 
are  the  proton  duration  x,and  variations  of  the  radio  frequency 
range  of  integration  such  as  oj2"  - or  the  full  available 

range.  It  is  not  clear  why  u^"-  does  not  work  as  well  as 

606-8800  MHz.  (See  Table  III). 

The  data  sample  of  just  19  events  cannot  be  considered  to 
be  large  enough  to  provide  decisive  answers  about  what  the  best 
parameters  or  procedures  are:  we  can  only  indicate,  from  this 
limited  data  base,  what  ought  to  be  explored  further.  In  this 
sense,  the  central  issue  we  have  raised  here  is  that  of  trying 
to  correlate  the  total  radio  energy  (for  the  U-time)  to  the  total 
proton  energy,  and  the  improved  correlations  suggest  that  this 
inquiry  should  be  pursued  further  over  an  extended  period  of  time. 
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